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Abstract. Two ultracompact Hll regions (IRAS 19181+1349 and 20178+4046) and one compact molecular clump 
(20286+4105) have been observed at far infrared wavelengths using the TIFR 1 m balloon-borne telescope and at 
mid infrared wavelengths using ISO. Far infrared observations have been made simultaneously in two bands with 
effective wavelengths of ~ 150 and ~210 /rm, using liquid ®He cooled bolometer arrays. ISO observations have 
been made in seven spectral bands using the ISOCAM instrument; four of these bands cover the emission from 
Polycyclic Aromatic Hydrocarbon (PAH) molecules. In addition, IRAS survey data for these sources in the four 
IRAS bands have been processed using the HIRES routine. In the high resolution mid infrared maps as well as 
far infrared maps multiple embedded energy sources have been resolved. There are structural similarities between 
the images in the mid infrared and the large scale maps in the far infrared bands, despite very different angular 
resolutions of the two. Dust temperature and optical depth (ri 5 o) maps have also been generated using the data 
from balloon-borne observations. Spectral energy distributions (SEDs) for these sources have been constructed 
by combining the data from all these observations. Radiation transfer calculations have been made to understand 
these SEDs. Parameters for the dust envelopes in these sources have been derived by fitting the observed SEDs. 

In particular, it has been found that radial density distribution for three sources is diffrent. Whereas in the case 
of IRAS 20178+4046, a steep distribution of the form r“^ is favoured, for IRAS 20286+4105 it is r“^ and for 
IRAS 19181+1349 it the uniform distribution (r'’). Line ratios for PAH bands have generally been found to be 
similar to those for other compact Hll regions but different from general Hll regions. 

Key words, infrared: interstellar: continuum - interstellar medium: H ll regions - interstellar medium: individual 
objects: (IRAS 19181+1349, IRAS 20178+4046, IRAS 20286+4105) 


1. Introduction 

Ultracompact Hii regions are a class of interesting as¬ 
tronomical sources that represent early stages of deeply 
embedded high mass (O or early B) stars. Therefore 
their study can provide vital information about high mass 
star formation as well as their interaction with the par¬ 
ent molecular cloud. They are also often associated with 
OH/H 2 O masers. Being deeply embedded in dust, almost 
all of their energy is emitted in the infrared waveband. 
In this paper we report far and mid infrared observations 
of two ultracompact Hii regions, IRAS 19181+1349 and 
IRAS 20178+4046. We also report observations for a com¬ 
pact molecular clump, IRAS 20286+4105, which is similar 
to the above sources in many respects. 

IRAS 19181+1349 and IRAS 20178+4046 are ultra¬ 
compact H II regions (Kurtz et al. 1994) with rising far in- 
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frared spectra between 60 and 100 /rm. Radio observations 
with the VLA exhibit multiple compact sources in the field 
of IRAS 19181+1349 (Kurtz et al. 1994; Zoonematkermani 
et al. 1990). Extended radio emission at 3.6 cm has also 
been seen in both the sources; however only for IRAS 
19181+1349 is the extended emission likely to be con¬ 
nected with the ultracomact Hii region (Kurtz et al. 
1999). Molecular lines of CO and CS have been observed 
towards both of these sources (Shepherd and Churchwell, 
1996; Bronfman et al. 1996). IRAS 19181+1349 is also as¬ 
sociated with H 2 O and OH masers (Forster and Caswell 
1989, Palagi et al. 1993); formaldehyde and methanol 
masers have not been detected (Caswell et al. 1995; 
Mehringer et al. 1995). No maser has been found to be 
associated with IRAS 20178+4046. Between 8 and 22 
^m IRAS low resolution spectrometer (LRS) spectrum is 
available for IRAS 20178+4046 (Volk and Cohen 1989) 
and groundbased intermediate resolution spectroscopic 
observations have been made between 3 to 13 /im (Faison 
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et al. 1998). In the latter, PAH features^ at 3.3 ^m, 8.7 
/rm and 11.3 /im, [Ne II] line at 12.8 /im and silicate fea¬ 
ture at 9.7 /rm have been detected. Based on radio re¬ 
combination observations, the kinematic distance of IRAS 
I9I8I-I-1349 has been variously determined as 9.7 kpc to 
11.8 kpc (Caswell et al. 1975, Kurtz et al. 1994). For IRAS 
20178-1-4046 the kinematic distances from 1.5 kpc to 3.3 
kpc have been used (Kurtz et al. 1994, Casoli et al. 1986). 
We use a distance of 9.7 kpc for IRAS 19181-1-1349 and 
3.3 kpc for IRAS 20178-f4046. 

IRAS 20286-1-4105 is an unusual compact CO source 
(Odenwald 1989) with infrared colours similar to ultra¬ 
compact H II regions (according to the criteria of Wood 
and Churchwell, 1989), associated with a young stellar 
object (YSO). It is located at the centre of a thin lu¬ 
minous ring of ~ ll5 diameter in the visible (Odenwald 
and Schwartz 1989). Other molecules like NH 3 , H 2 O and 
CS have also been observed towards this source (Molinari 
et al. 1996, Bronfman et al. 1996). Near infrared imag¬ 
ing of this source in J, H and K bands has been done 
by Comeron and Torra (2001). LRS spectrum of IRAS 
20286-1-4105 (Olnon et al. 1986) shows a broad absorption 
feature at 11.3 /im. It has also been associated with H 2 O 
masers (Codella et al. 1996). The kinematic distance for 
this source is 3.7 kpc (Molinari et al. 1996). 

As part of a programme to obtain far infrared images 
of star forming regions (Ghosh et al. 1988, 1989a, 1989b, 
1990, 2000; Mookerjea et al. 1999) we have mapped these 
three regions simultaneously in two far infrared bands us¬ 
ing the 1 m balloon-borne telescope of the Tata Institute of 
Fundamental Research (TIFR). The effective wavelengths 
of the two bands are ~ 150 /im and ^ 210 /im and the 
angular resolution acheived is ^ V at both wavelengths. 
We have also procured HIRES (Aumann et al. 1990) pro¬ 
cessed IRAS maps for these sources in all the four IRAS 
bands. To understand the role of transient heating of PAH 
molecules, we have imaged the central 3' x 3' regions of 
these sources in seven mid-infrared bands using ISOCAM 
in the ISO| satellite. 

The observational details are given in Sect. 2. The re¬ 
sults are presented in Sect. 3. Spectral energy distribu¬ 
tions (SEDs) have been constructed by combining our 
results with other observations. To understand the in¬ 
frared emission from these sources, radiation transfer cal¬ 
culations have been performed. The schemes of radiation 
transfer calculations are presented in Sect. 4. Discussion 
of sources is presented in Sect. 5. Conclusions are summa¬ 
rized in Sect. 6 . Preliminary results of these observations 
have been given in Verma et al. (1999a,b) and Karnik 
( 2000 ). 


^ These features are also called unidentified infrared bands 
(UIBs) meaning that the PAH origin is not completely settled. 
However we shall consistently use the term PAH. 

^ Based on observations with ISO, an ESA project with in¬ 
struments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United 
Kingdom) with the participation of ISAS and NASA. 


2. Observations 

2.1. Balloon-borne observations 

Balloon-borne far infrared observations were made with 
the TIFR I m telescope during a balloon flight on 
February 20, 1994 from Hyderabad (lat. 17?47 N, long. 
78?57 E), India. The details of the telescope and the ob¬ 
servational procedures have been described by Ghosh et 
al. (1988); here only some basic features will be mentioned. 
The telescope is a 100 cm aperture f /8 Cassegrain system. 
It can be pointed towards a desired direction in the sky 
with an rms pointing stability of better than ~ 0f3. The 
incident radiation is chopped by vibrating the secondary 
mirror along the cross elevation direction at a frequency 
of 10 Hz and a throw of 4.'2. The detectors consist of two 
arrays (each 2x3) of composite Si bolometers cooled to 
0.3 K by a closed cycle ^He refrigerator. The field of view 
for each detector is 1.'6. The observations were made si¬ 
multaneously in two bands with effective wavelengths of 
^ 150 /im and 210 /im. The two bands are separated by 
the use of cool Csl restrahlen beam-splitter which reflects 
the radiation with wavelengths between 85 /tm and 180 /tm 
and transmits the radiation with wavelengths shorter than 
85 /im and longer than 180 /im. The radiation with wave¬ 
lengths shorter than 85 /tm is blocked by other blocking 
filters. Some of the details of the photometer are given in 
Verma et al. (1993). Measured transmissions for both the 
bands as well as the effective wavelengths and the relative 
responsivities (which depend on the source spectrum) are 
given in Ghosh et al. (2000). For a dust temperature of 
30 K and a gray body spectrum with emissivity, the 
effective wavelenghts for the two bands are 148 /tm and 
209 /im. For simplicity, we convert the flux densities in 
the two bands to values at wavelengths of 150 and 210 
/tm respectively. 

Mapping was done in a raster mode with a scan speed 
of 0f8 along the cross elevation axis. The boresight of 
the telescope was determined for each signal sample using 
the data from the star tracker and the rate gyroscopes. For 
each band, signals from all six detectors were used to form 
a two dimensional signal grid in instrument coordinates 
(elevation, cross-elevation) with a cell size of 0f3 x 0f3. 
Data have been deconvolved using a Maximum Entropy 
Method (MEM) similar to that of Gull and Daniell (1978). 
Absolute positions were calibrated using observations of 
catalogued stars with an optical photometer in the focal 
plane; field of which is offset with respect to the infrared 
field. The calibrated absolute positions are accurate to 
V. Jupiter was observed at the begining of the flight as 
well as towards the end of the flight; Jupiter observations 
were used for flux calibration. As the Jupiter size (0.'6) was 
much less than the field of view (116), Jupiter observations 
were also used for the determination of the point spread 
function (PSF) for MEM deconvolution. For the purpose 
of flux calibration we have assumed Jupiter to be a black 
body with a temperature of 127 K. The final resolution 
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in the deconvolved maps as given by the FWHM (El. x 
XEl.) of Jupiter is l.'O x 1.'3 in both the bands. 

Taking advantage of the simutaneous observations in 
the two bands, with almost identical field of view, we have 
generated maps of the dust temperature (T^;) and optical 
depth at 150 /im (ri5o). The procedure for the generation 
of these maps is described in Mookerjea et al. (2000). For 
these maps we have assumed the dust emissitivity of ex oc 
and a binning of 3 pixel by 3 pixel has been effected 
for both the bands before computing the Td and nso. The 
temperature and optical depth maps have been restricted 
to regions where the flux densities in both the maps are 
well within the dynamic range. 

2.2. IRAS observations 

IRAS survey data for these three sources were processed 
using the HIRES routine at the Infrared Processing and 
Analysis Center (IPAC), Caltech ^ to obtain high res¬ 
olution images in all the four IRAS bands. Details of 
the HIRES routine have been given in Aumann et al. 
(1990). Angular resoluions in different images along ob¬ 
served source sizes are given in Table 1. The HIRES maps 
have been used to obtain flux densities and angular sizes 
in the IRAS bands. 

2.3. ISO observations 

The central portions of the sources were observed at mid 
infrared wavelengths using ISOCAM. An area of 3' x 3' 
around each source has been imaged in seven spectral fil¬ 
ters. The filters are SW2, (with central wavelengths of 3.30 
/rm; bandwidths of 0.20 ^m), SW6 (3.72; 0.55), LW4 (6.00; 
1.0), LW5 (6.75; 0.5), LW6 (7.75; 1.5), LW7 (9.62; 2.2) and 
LW8 (11.4; 1.3) (ISOCAM Observer’s Manual, 1994). Out 
of these filters, four (SW2, LW4, LW6 and LW8) include 
emission from PAH bands whereas the remaining three 
bands are used as comparison for continuum emission. The 
pixel size used for imaging was &' x 6". Integration time 
per exposure was 2.0 s for SW filters and 0.28 s for LW 
filters. A sufficient number of exposures were taken for 
detector stabilization before the actual observations. The 
data were analyzed using Cam Interactive Analysis (CIA) 
version 3.0. 

3. Results 

3.1. IRAS 19181+1349 

Around IRAS 19181-1-1349 an area of ~ 24' x 30' was 
scanned twice and the resulting deconvolved maps in the 
two far infrared bands are shown in Fig. 1. For this re¬ 
gion no suitable star was seen by the optical photome¬ 
ter, therefore the expected error in absolute position is 
larger (~ 2'). We have made corrections to the coordi¬ 
nates to match the position of the main peak with that 

® IPAC is funded by NASA as part of the IRAS extended 
mission program nnder contract to JPL 


in the IRAS HIRES map at 100 /rm. As can be seen from 
the figure, there are multiple sources in both the maps. 
Whereas in the longer wavelength band the source is re¬ 
solved into three components and a lobe in the west, in 
the shorter wavelength band it is resolved into two compo¬ 
nents. The deconvolved FWHM size of IRAS 19181-1-1349 
at 150 fim is 2.'1 x 317 and at 210 /rm it is 212 x 213 ; the lat¬ 
ter encloses the three components. The second component 
seen at 150 fim is not covered by ISOCAM observations. 
Besides IRAS 19181-I-1349, many other sources are seen 
in the maps. Among these are sources corresponding to 
IRAS 19175-bl357 and IRAS 19175-fl342. Also, there is 
a lot of diffuse emission in both the bands. Integrated flux 
densities around the peaks are given in Table 2. 

The dust temperature (T^;) and optical depth (tisq) 
maps for IRAS 19181-1-1349 are shown in Fig. 2. It is seen 
that these maps are complex. There are multiple peaks in 
temperature which are shifted with respect to the main 
intensity peak by 1' to 4'. Optical depth map has three 
peaks. 

The HIRES processed maps of IRAS 19181-I-I349 in 
the four IRAS bands, derived from the survey data are 
shown in Fig. 3. In HIRES maps also the source is resolved 
at all four wavelengths. It shows two components at 12 and 
25 /im but at 60 and 100 /tm only one component is seen. 
However, as seen from Table 1, even at 60 and 100 /im the 
source is extended with FWHM of 1.'3 x 2H and 213 x 311 
respectively as compared to the resolution of 111 x 116 and 
210 X 213. 

Maps generated from ISOCAM images in the seven 
bands are shown in Fig. 4. The images in the SW filters 
have been physically shifted such that the positions of the 
main peaks in these images coincide with the respective 
positions in the LW filters. This was done to correct for 
the lens wheel jitter. This procedure has also been followed 
for the SW images of the other two sources. The required 
shift is ~ 20". It can be seen from this figure that IRAS 
19181-1-1349 consists of multiple components in all the fil¬ 
ters. However the secondary peaks in the two SW filters 
do not appear in the LW filters and therefore do not seem 
to be associated with the IRAS source. The components 
in the LW filters can be divided into two groups separated 
by ^ 115, mainly in the E-W direction. This feature is 
also seen at 12 and 25 /im in the HIRES maps and at 210 
/im in the TIER map. Integrated flux densities of strong 
components are given in Table 3. 

3.2. IRAS 20178+4046 

Around IRAS 20178-1-4046 an area of ~ 30' x 30' was 
scanned twice and the resulting deconvolved maps in the 
two far infrared bands are shown in Fig. 5. This source 
does not show multiple components in any of the two 
maps. It can be seen from Table 1 that the source is partly 
resolved. 

The dust temperature (T^) and optical depth (tisq) 
maps for IRAS 20178-1-4046 are shown in Fig. 6. It is seen 
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Table 1. Map resolutions and source sizes 


Source 


FWHM Sizes in arcminutes (minor x major) 


HIRES-IRAS 


TIFR 

12 /xm 

25 /xm 

60 /xm 

100 /xm 

150 /xm 

210 /xm 

19181+1349 

Map resolution 

0.5x1.2 

0.5X1.0 

1.1x1.6 

2.0x2.3 

+0x1.3 

+0x1.3 


Source size 

+3x2.8 

0.9x2.1 

1.3x2.1 

2.3x3.! 

2.1x3.7 

2.2x2.3 

20178+4046 

Map resolution 

0.5x1.1 

0.6x1.2 

1.1x1.8 

1.9x2.4 

+0x1.3 

+0x1.3 


Source size 

0.6X0.6 

0.5X0.7 

1.0x2.0 

1.9x2.2 

+1x1.5 

+5x1.5 


Model size 

0.6X1.3 

0.6X0.6 

1.1x2.! 

1.9x2.5 

+2x1.4 

+2x1.5 

20286+4105 

Map resolution 

0.5X1.0 

0.6x1.2 

1.0x1.1 

1.9x2.4 

+0x1.3 

+0x1.3 


Source size 

0.8X1.3 

0.6X0.9 

1.1x1.5 

1.8x2.0 

0.9x1.2 

+2x1.5 



Fig. 1. Intensity maps of the region around IRAS 19181+1349 at (a) 150 /im and (b) 210 /rm. The dotted boundary 
marks the area scanned. Crosses mark the positions of IRAS PSC sources. The contour levels are at 0.90, 0.70, 0.50, 
0.30, 0.20, 0.10, 0.05 and 0.025 of the peak intensity. The peak intensities are 723 Jy arcmin“^ at 150 /im and 390 Jy 
arcmin“^ at 210 /xm. 


that the peak in optical depth occurs close to the peak 
in the intensity distributions. The temperature distribu¬ 
tion shows two peaks separated by 2!4. The total range in 
temperature is from 72 K to 21 K. 

The HIRES processed maps of IRAS 20178+4046 in 
the four IRAS bands, derived from the survey data, are 
shown in Fig. 7. The central part of the source is not 
resolved in any of the IRAS bands, as seen from Table 1. 
However there is lower level extended emission in all the 
bands except at 25 /im. 

Maps generated from the ISOCAM images in the seven 
bands are shown in Fig. 8. It can be seen from this figure 


that IRAS 20178+4046 consists of only one component 
in all the filters. The contours in all the LW filters are 
remarkably similar to each other, showing a lobe towards 
the south. 

3.3. IRAS 20286+4105 

Around IRAS 20286+4105 an area of ~ 16' x 30' was 
scanned twice and the resulting deconvolved maps in the 
two far infrared bands are shown in Fig. 9. As in the case 
of IRAS 20178+4046, this source also has only one com¬ 
ponent in both the bands. The source is resolved at 210 
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Table 2. Far infrared flux densities of the sources 


Source RA Dec. Flux Density in Jy* 

HIRES-IRAS TiFR 



(1950) 

(1950) 

12 /im 

25 /tm 

60 /tm 

100 /im 

150 /im 

210 /im 

19181+1349 

19 18 09.6 

13 49 46 

156 

669 

6443 

7770 

3939 

2590 

19181+1349E 

19 18 12.7 

13 49 55 

46 

267 

<2427 

<2236 

<1096 

656 

19181+1349W 

19 18 07.4 

13 49 39 

49 

265 

<1796 

<1679 

<1496 

827 

19176+1359 

19 17 36.8 

13 59 16 

14 

47 

<473 

<846 

572 

351 

19175+1357 

19 17 32.7 

13 57 16 

22 

95 

<877 

1254 

<680 

<367 

20286+4105 

20 28 38.5 

41 05 43 

39 

170 

1071 

1348 

1249 

902 

20178+4046 

20 17 48.6 

40 46 44 

114 

748 

2932 

3325 

2244 

1199 


*The flux densities have been integrated in a circle of 4' diameter except for IRAS 19181+1349 where integration has been 
done in a circle of 5' diameter; for the two components 19181+1349E and 19181+1349W the integrations have been done in a 
circle of 2' diameter. 



Fig. 2. Dust temperature (T^) and optical depth (nso) distribution for the region around IRAS 19181+1349. (a) 
Temperature distribution. The contours are at the levels of 50, 40, 35, 30, 25, 20 and 15 K. (b) Optical depth 
distribution. The contours are at 90, 70, 50, 30, 20, 10 and 5% of the peak value which is 0.1. 


/im but at 150 fxm it is unresolved. At 210 /rm it shows an 
extension in the same direction (N-S) in which the source 
shows double structure in ISO maps (see later). 

The dust temperature (T^) and optical depth (tiso) 
maps for IRAS 20286+4105 are shown in Fig. 10. The dust 
in IRAS 20286+4105 is very much cooler (temperature 
range of 15 K to 25 K) than the other two sources. This 
is perhaps because the luminosity for this source is lowest 
(see Table 4). The peak in temperature occurs near the 
peak of the flux density distribution. However the peak in 
optical depth is shifted to the north by ^ 1'. 

The HIRES processed maps of IRAS 20286+4105 in 
the four IRAS bands, are shown in Fig. 11. As seen from 
Table 1, the source is only resolved at 12 /xm and partly 
resolved at 60 /im. 


Maps generated from the ISOCAM images in the seven 
bands are shown in Fig. 12. It can be seen from this figure 
that IRAS 20286+4105 consists of two main components, 
separated along the N-S direction by ~ 32" in all the 
filters. There is a third weaker component to the east. The 
structure in our maps (especially the SW maps) is very 
similar to the structure in the K band image of Comeron 
and Torra (2001). 

4. Radiation Transfer Calculations 

To understand the infrared emission from these sources 
radiation transfer calculations have been performed. The 
details of the basic radiation transfer scheme are given in 
Mookerjea et al. (1999). The geometry of the basic model 
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Fig. 3. The HIRES processed intensity maps of the region around IRAS I9I8I+1349 at (a) 12 /im, (b) 25 /rm, (c) 
60 /rm and (d) 100 ^m. The contour levels are at 0.90, 0.70, 0.50, 0.30, 0.20, 0.10, and 0.05 of the peak intensity for 
12 /im; for other wavelengths, besides these there are additional contours at 0.025 of the peak. The peak intensities 
are 26.4, 181, 1580 and 1090 Jy arcmin“^ at 12, 25, 60 and 100 /tm respectively. 


is shown in Fig. 13a. Briefly, an internal radiation source, 
i.e., a star is surrounded by a spherically symmetric en¬ 
velope of dust and gas mixture. The gas is assumed to 
be hydrogen and surrounds the star from the stellar sur¬ 


face itself. The dust is a mixture of astronomical silicates 
and graphite and there is an inner cavity of radius Rmm 
which is devoid of dust. Two types of interstellar dust are 
considered; those with absorption and emission properties 
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Table 3. Flux densities of sources in the ISOCAM images 


IRAS 

Source 

RA 

(1950) 

Dec. 

(1950) 



Flux Density in 

Jy* 



SW2 
3.30 /im 
PAH 

SW6 
3.72 /xm 

LW4 
6.00 /xm 
PAH 

LW5 
6.75 /xm 

LW6 
7.75 /xm 
PAH 

LW7 
9.62 /xm 

LW8 
11.4 /xm 
PAH 

19181-kl349 

19 18 05.4 

+13 49 22 

- 

- 

0.81 

1.34 

2.6 

- 

1.29 


19 18 05.8 

+13 50 32 

O.IO’^ 

- 

- 

- 

- 

- 

- 


19 18 06.8 

+13 49 40 

- 

- 

0.73 

1.22 

2.4 

1.29 

3.0 


19 18 07.6 

+13 48 54 

0.17# 

0.19# 

- 

- 

- 

- 

- 


19 18 11.2 

+13 49 33 

0.19 

- 

1.21 

1.70 

3.6 

- 

1.5 


19 18 12.2 

+13 50 52 

0.09 

- 

- 

- 

- 

- 

- 


19 18 12.5 

+13 51 24 

0.21 

0.20 

- 

- 

- 

- 

- 


19 18 12.6 

+13 50 10 

0.77 

0.56 

2.5 

2.8 

5.9 

1.90 

3.4 


19 18 13.0 

+13 49 59 

0.23# 

0.13# 

1.27# 

1.51# 

2.85# 

- 

- 


19 18 13.5 

+13 49 24 

- 

- 

0.78 

0.96 

2.0 

- 

- 

20178-74046 

20 17 53.3 

+40 47 03 

1.8 

0.70 

11.6 

10.7 

30.4 

11.0 

20.9 


20 17 55.0 

+40 47 01 

0.30# 

- 

2.7# 

2.7# 

7.5# 

2.7# 

4.6# 

20286+4105 

20 28 40.3 

+41 05 14 

0.71 

0.23 

4.4 

4.8 

11.5 

3.5 

6.3 


20 28 40.3 

+41 05 45 

0.25 

0.13 

3.1 

3.3 

7.8 

1.90 

2.8 


20 28 42.5 

+41 05 44 

0.20 

0.12 

2.9 

2.4 

3.5 

1.21 

1.50 


*The flux densities have been integrated in a circle of 30” diameter around the peak except for those marked with where the 
integration is done over an area of 18” x 18” around the peak. 


given by Draine & Lee (1984) (hereafter DL) and those 
with properties given by Mathis et al. (1983) (hereafter 
MMP) and the size distribution of the grains is taken to 
be a power law given by Mathis et al. (1977), viz. n(a)da 
(X a“™da, with m = 3.5 and 0.01 /xm < a < 0.25 /xm; a 
is the radius of the grain. The envelope is surrounded by 
the interstellar radiation field (ISRF) taken from Mathis 
et al. (1983). The radial density distribution of the gas 
and dust has been taken to be a power law, i.e. n(r) oc 
r““ with a = 0, 1 or 2. The radiation transport has been 
carried out using a program based on the code CSDUST3 
(Egan et al. 1988). The parameters of the model are - 
inner radius of the dust cloud Rmi„, outer radius of the 
envelope Rmax, a, luminosity and temperature of the em¬ 
bedded star, relative abundance of the two types of dust, 
gas to dust ratio and the total radial optical depth at a 
selected wavelength (noo at 100 ^m). The radiation trans¬ 
fer of Lyman continuum photons in gas is carried out in a 
consistent manner with that through the dust. The com¬ 
bination of the parameters is obtained which gives lowest 

to fit the observed SED. The scheme also gives the 
expected radio emission and the radial distribution of in¬ 
frared emission. These have also been used to constraint 
the model. This basic scheme was used to model IRAS 
20178-k4046. 

As ISOCAM maps of IRAS 20286-1-4105 show two 
main components in all the images, the radiation trans¬ 
fer calculations for this source were performed differently. 
Whereas the basic scheme of the radiation transfer is the 
same, it is assumed that there are two cores. Each core, 
with its own internal radiation source, is surrounded by a 


spherically symmetric envelope of dust and gas mixture. 
The type of dust (DL or MMP) and the silicate to graphite 
ratio is assumed to be the same for both cores. Total lu¬ 
minosity of the source is divided between the two cores in 
the ratio of their respective ISOCAM flux densities. The 
radiation transfer calculation is done for each core sepa¬ 
rately and the radiation emerging from each core is cal¬ 
culated. These cores are surrounded by another common 
shell, with the internal input radiation equal to the sum 
of the emerging radiation from each core. This common 
shell is surrounded by interstellar radiation field taken 
from Mathis et al. (1983). The radiation emerging from 
the inner boundary of this shell is the input at the outer 
boundary of the cores. The calculations are done itera¬ 
tively. The geometry of the model is shown in Fig. 13b. 

The radiation transfer calculations for IRAS 
19181-1-1349 were performed in a cylindrical geome¬ 
try with two embedded sources. This model has been 
described in Karnik and Ghosh (1999). An attempt was 
made to fit the SED as well as the radial profile of the 
intensities. 

5. Discussion 

5.1. IRAS 20178+4046 

First we consider IRAS 20178-1-4046 which is the simplest 
of the three sources because it shows only one component 
in all the wavelength bands. We have constructed the SED 
of IRAS 20178-1-4046 using the flux densities from balloon- 
borne, IRAS and ISO observations, given in Tables 2 and 
3, IRAS LRS data from Volks and Cohen (1989) and 
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Fig. 4. The maps of ISOCAM observations of the region around IRAS 19181+1349 at (a) 3.30 /rm, (b) 3.72 /xm, 
(c) 6.00 ^m, (d) 6.75 ^m, (e) 7.75 ^m, (f) 9.62 ^m, and (g) 11.4 ^m. The contour levels are at 0.90, 0.70, 0.50, 
0.30, 0.20, 0.10 and 0.05 of the peak intensity. The peak intensities are 19.5, 17.8, 36.5, 37.7, 70.2, 37.8 and 58.7 Jy 
arcmin”^ at 3.30, 3.72, 6.00, 6.75, 7.75, 9.62 and 11.4 respectively. 


ground based observations of Faison et al. (1998). This 
SED is shown in Fig. 14a. In this figure, we have only 
shown ISOCAM flux densities in the continuum bands; 
PAH bands are treated separately. It can be seen that 
while there is general agreement between various observa¬ 
tions, HIRES flux density at 12 /xm is higher than IRAS 
LRS data as well as the data of Faison et al. (1998). This is 
because the HIRES flux density is integrated over 4' dia 
and includes substantial amount of extended flux. Also, 
the flux densities of Faison et al. (1998) are lower than 
those of IRAS LRS as well as from ISO observations. This 
is because of the small aperture (9") of Faison et al. (1998). 

The SED from the best fit radiation transfer model is 
shown in Fig. 14a. It can be seen that the calculated SED 


fits the observed one quite well. The parameters of the 
best fit model are given in Table 4. The luminosity cor¬ 
responds to an 07 ZAMS star. The density distribution 
that fits best is, n(r) oc r“^, the dust type is MMP and is 
dominated by silicates. The calculated angular sizes (de¬ 
convolved with the angular resolutions) are given in Table 
1. It can be seen that the calculated angular sizes are 
consistent with the observed sizes. The calculated radio 
flux density at 6 cm is 72 mJy for a gas to dust ratio of 
100. This can be compared with the observed flux density 
of 69 mJy by McCutcheon et al. (1991) and 65 mJy by 
Wilking et al. (1989). Density distribution (r“^) is con¬ 
sistent with what is expected in the envelopes of young 
protostars from free fall models (see for example Shu et 
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Fig. 5. Same as figure 1 but for the region around IRAS 20178+4046. The peak intensities are 978 Jy arcmin ^ at 
150 /xm and 589 Jy arcmin”^ at 210 /xm. 



Fig. 6. Same as figure 2 but for the region around IRAS 20178+4046. The temperature contours are at the levels of 
50, 40, 35, 30 and 25 K and for the optical depth distribution the peak value is 0.038. 


al. 1987). While it is true that the fitting of SED alone is 
not unique (Churchwell et al. 1990), in the present case 
angular distribution and radio flux are consistently mod¬ 
elled and therefore parameters of the model are on firmer 
footing. 


5.2. IRAS 20286+4105 

We have constructed the SED for IRAS 20286+4105 us¬ 
ing the flux densities from balloon-borne, IRAS and ISO 
observations, given in Tables 2 and 3 and IRAS LRS data 
from Olnon et al. (1986). This SED is shown in Fig. 14b. 
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Fig. 7. The HIRES processed intensity maps of the region around IRAS 20178+4046 at (a) 12 /rm, (b) 25 /im, (c) 60 
/rm and (d) 100 ^m. The contour levels are at 0.90, 0.70, 0.50, 0.30, 0.20, 0.10, 0.05 and 0.025 of the peak intensity; 
for 25 /im additional contour at 0.01 of the peak has been shown. The peak intensities are 89.7, 1490, 1110 and 737 
Jy arcmin”^ at 12, 25, 60 and 100 /im respectively. 


For ISO observations we have combined the flux densities 
of two components (N and S); the third component toward 
the East gets weaker at longer wavelengths and therefore 
has not been included. It can be seen that there is good 
agreement between various observations. 

The SED from the best fit radiation transfer model is 
shown in Fig. 14b. The calculated SED fits the observed 
SED quite well. The parameters of the best fit model are 
given in Table 4. Here the two cores are denoted as IRAS 
20286+4105 N and IRAS 20286+4105 S. The luminosi¬ 
ties of the two cores correspond to ZAMS stars of the 
type BO.5 (N) and BO-0.5 (S). The dust density distri¬ 
bution that fits best is, n(r) oc r“^ and the dust type 
is again MMP and dominated by silicates. The radiation 
transfer model predicts negligible radio emission. No com¬ 


pact radio source has been observed in the field of IRAS 
20286+4105. McCutcheon et al. (1991) see an extended 
feature at 6 cm and according to them the association 
of radio continuum with IRAS 20286+4105 is not clear. 
Similarly Wendker et al. (1991) see radio emission from 
IRAS 20286+4105 as a small density enhancement on a 
large ridge. Thus the model is consistent with radio obser¬ 
vations. The flatter dust density distribution, as compared 
to IRAS 20178+4046, suggests that IRAS 20286+4105 is 
relatively more evolved system. 

5.3. IRAS 19181+1349 

The SED for IRAS 19181+1349 using the flux densities 
from balloon-borne, IRAS and ISO observations, is shown 
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Fig. 8. Same as figure 4 but for the region around IRAS 20178+4046. The peak intensities are 30.3, 15.2, 232, 214, 
617, 239 and 505 Jy arcmin”^ at 3.30, 3.72, 6.00, 6.75, 7.75, 9.62 and 11.4 /xm respectively. 


Table 4. Parameters* for the best fit radiation transfer models 


Source 

Rmacc 

(pc) 

Rmin 

(pc) 

L 

(10" L«) 

a 

noo 

Dust Type 

^Silicate 

20178+4046 

3.40 

0.077 

9.5 

2 

0.12 

MMP 

0.76 

20286+4105 N 

0.24 

0.022 

1.2 

1 

0.11 

MMP 

0.65 

20286+4105 S 

0.32 

0.022 

2.3 

1 

0.11 

MMP 

0.65 

20286+4105 Outer 

3.02 

1.12 

3.5 

1 

0.025 

MMP 

0.65 

19181+1349 

1.8 

0.03, 0.01 

62, 63 

0 

0.10 

DL 

0.80 


*The scheme of calculation for IRAS 19181+1349 is different from that of other sources. Therefore some of the parameters are 
not directly comparable. 
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Fig. 9. Same as figure 1 but for the region around IRAS 20286+4105. The peak intensities are 770 Jy arcmin ^ at 
150 /xm and 402 Jy arcmin”^ at 210 /xm. 


in Fig. 14c. As was mentioned earlier, maps of IRAS 
19181+1349 show multiple components in all the images. 
For this source we have used the flux density within 5^ of 
the main peak such that the contribution from multiple 
sources is included. 

The SED for IRAS 19181+1349 was analyzed by 
Karnik and Ghosh (1999) using a model with two em¬ 
bedded sources in a cylindrical geometry. An attempt was 
made to fit the SED as well as radial profile of the in¬ 
tensities. The best fit to the SED is shown in Fig. 14c. 
The parameters of that model are given in Table 4. The 
prefered dust density distribution (r°) suggests that IRAS 
19181+1349 is very evolved system. 

5.4. PAH emission 

As seen from Table 3, the emission in PAH bands is gen¬ 
erally higher as compared to the neighbouring bands. We 
have derived line fluxes in the PAH bands from the flux 
densities given in Table 3. For this, we have subtracted 
the continuum flux density taken from the neighbouring 
filter; SW6 (3.72 /xm) for 3.30 ^m band, LW5 (6.75 ^m) 
for 6.2 and 7.7 /xm bands and LW7 (9.6 /xm) for 11.3 /xm 
band; and mutiplied by the filter bandwidth. This proce¬ 
dure will lead to an overestimation for the 11.3 ^m PAH 
feature because the continuum given by LW7 is underes¬ 
timated due to silicate absorption. The band fluxes and 


the band ratios are given in Table 5. It is seen that the 
strongest PAH feature is the one at 7.7 /xm whereas the 
feature at 6.2 /xm is absent in most of the cases. We have 
also given the range of band ratios obtained by Roelfsema 
et al. (1996) using ISO Short Wavelength Spectrometer 
(SWS) for six compact H ii and the average value of line 
ratios from typical Hii regions from Cohen et al. (1989). 
It is seen that the band ratios obtained by us are consis¬ 
tent with the values of Roelfsema et al. (1996) for three 
of the four features observed by us. However they are dif¬ 
ferent from the values of Cohen et al. (1989) for typical 
Hii regions. For example, the ratio of 7.7 /xm and 3.3 /xm 
features varies from 18 to 52 in our data as compared to 
17 to 73 in those of Roelfsema et al. (1996), but the av¬ 
erage value of Cohen et al. (1989) is 7. This shows that 
there is a difference in the PAH emission by compact H ii 
regions and other H ii regions. However, we can not un¬ 
derstand the near absence of the 6.2 /xm PAH band in our 
observations especially as both, the 6.2 /xm feature and 
the strongest observed feature at 7.7 /xm arise from C-C 
stretching and have been found to be strongly correlated 
(Cohen et al. 1989). 

In addition, an attempt has been made to generate 
the spatial distribution of emission in individual PAH fea¬ 
tures for the regions imaged by ISOCAM filters. The PAH 
emission maps have been generated by subtracting out the 
expected continuum emission. The spectrally local contin- 
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Fig. 10. Same as figure 2 but for the region around IRAS 20286+4105. The temperature contours are at the levels of 
25, 22, 20, 18 and 16 K and for the optical depth distribution the peak value is 0.042. 


uum emission has been estimated by linear interpolation 
or extrapolation using the neighbouring ISOCAM filters 
(viz., LW5 & LW7 for the features at 7.7 /rm and 11.3 /xm 
features; LW5 & SW6 for the features at 3.3 ^m and 6.2 
fxm). Generally the above method has been successful, ex¬ 
cept for the case of 3.3 ^m feature for IRAS 20178+4046 
and IRAS 20286+4105. The failure of the latter cases is 
attributed to a large ratio (~ 15) of intensities between 
6.75 /im and 3.72 /xm, making linear extrapolation a sus¬ 
pect. In these cases the measured intensities in the nearest 
filter (SW6) itself have been used as an estimate of the lo¬ 
cal continuum. 

The spatial distribution of emission in individual PAH 
features is found to be diffuse and widespread in na¬ 
ture. However, they qualitatively (morphologically) re¬ 
semble the isophotes of total emission in the correspond¬ 
ing ISOCAM filters (i.e. uncorrected for the continuum), 
which have been presented in figures 4, 8 and 12. Hence 
these are not presented here. Instead, some quantita¬ 
tive results are presented. The total emission in the in¬ 
dividual PAH features from the imaged regions of IRAS 
19181+1349, 20178+4046 and 20286+4105 are tabulated 
in Table 6. 

Considering that the ratios of emission in the individ¬ 
ual PAH features may contain information about their ex¬ 
citation mechanisms, these have been computed on pixel- 
by-pixel basis for the three sources considered here. The 
determination of these ratios have been restricted to re¬ 
gions with sufficient signal strengths in relevant PAH fea¬ 
tures. The summary of these ratios are presented in Table 
7. 


Table 7. Distribution of PAH feature ratios for the ISO 
images 


IRAS 

Source 

hi.sIh.T 

l7.7/h.2 

77 . 7 / 73.3 

19181+1349 

0.30+0.39* 

- 

23.2+11.8* 


(1024)“ 


(476)“ 

20178+4046 

0.84+0.40 

14.7+8.3 

19.9+14.0 


(289) 

(199) 

(143) 

20286+4105 

1.07+0.64 

- 

13.4+11.7 


(256) 


(145) 


‘mean ± std. dev. 

“(sample size : # of pixels) 


The scatters in the values of the feature ratios (see the 
standard deviations) is not measurement noise but repre¬ 
sent their genuine spatial variation. This implies different 
physical conditions around each of the star forming re¬ 
gions which differentially affect the exciting mechanisms 
of individual PAH features. Whereas the mean ratios are 
similar for the regions IRAS 20178+4046 and 20286+4105, 
the values for the region IRAS 19181+1349 differ signifi¬ 
cantly. This may be related to their evolutionary stages, 
which is also indicated by the conclusions from the ra¬ 
diative transfer modelling (see Sect. 4.1.3; viz., the power 
law index of radial density distribution is flattest for IRAS 
19181+1349, indicating it to be the most evolved source). 
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Table 5. PAH band fluxes for the sources in ISOCAM images 


IRAS 


RA 

Dec. 

Band Flux“ (10“ 

erg cm’ 

- s-“) 

PAH Band Ratios 

Source 

(1950) 

(1950) 

3.3 /rm 

6.2 

7.7 /rm 

11.3 fim 

6.2/3.3 

7.7/3.3 

11.3/3.3 

19181-tl349 

19 

18 05.4 

+13 49 22 

- 

- 

9.8 

- 

- 

- 

- 


19 

18 06.8 

+13 49 40 

- 

- 

9.1 

5.2 

- 

- 

- 


19 

18 11.2 

+13 49 33 

- 

- 

14.4 

- 

- 

- 

- 


19 

18 12.5 

+13 51 24 

0.06 

- 

- 

- 

- 

- 

- 


19 

18 12.6 

+13 50 10 

1.16 

- 

23 

4.6 

- 

20 

3.9 


19 

18 13.0 

+13 49 59 

0.55 

- 

10.1 

- 

- 

18 

- 


19 

18 13.5 

+13 49 24 

- 

- 

7.9 

- 

- 

- 

- 

20178-t4046 

20 

17 53.3 

-t40 47 03 

6.1 

7.6 

149 

30 

1.25 

25 

5.0 


20 

17 55.0 

-t40 47 01 

- 

- 

36 

5.8 

- 

- 

- 

20286-t4105 

20 

28 40.3 

+41 05 14 

2.6 

— 

51 

8.5 

— 

19 

3.2 


20 

28 40.3 

+41 05 45 

0.66 

- 

34 

2.6 

- 

52 

4.0 


20 

28 42.5 

+41 05 44 

0.44 

4.2 

8.3 

0.91 

9.5 

19 

2.1 

Compact H ll 

h 








3.4 - 22 

7-44 

2.9 - 6.7 


Typical H ll 

4.0 

7 

2.3 

regions“ 





“Band fluxes have been caculated after subtracting the continuum defined by the neighbouring filter as given in the text 
*'The flux ratios for compact Hll regions from Roelfsema et al. (1996) 

“The flux ratios for typical Hll regions from Cohen et al. (1989) 


Table 6. Total emission in individual PAH features extracted from ISOCAM data 


PAH 

feature 

Total 

emission in the PAH feature 
erg sec“^ cm“^ 


IRAS 19181-11349 

IRAS 20178-14046 

IRAS 20286-14105 

3.3 /im 

1.56 xlO"^ 

1.15 xl0“^“ 

5.20 xlO""^ 

6.2 /im 

- 

3.44 xl0“^° 

- 

7.7 /im 

6.03 xl0“® 

4.01 xlO"® 

1.56 xl0“® 

11.3 /tm 

7.45 xl0~“ 

9.62 xl0“^° 

3.60 xl0~“ 


6. Conclusions 

We have mapped two ultracompact H ii regions and one 
molecular clump simultaneously in two far infrared bands 
using 1 m TIFR balloon-borne telescope. Using the maps 
from the balloon-borne observations, maps of dust tem¬ 
perature and optical depth have been obtained. We have 
also imaged the central 3' x 3' regions of these sources 
in seven mid infrared bands using ISOCAM instrument of 
ISO. IRAS HIRES processed data have also been obtained 
for these sources in the four IRAS bands. There are certain 
similarities between various maps of the sources. Using the 
flux densities in all the 13 bands, as well as any other avail¬ 
able data, SEDs have been constructed. Radiation trans¬ 
fer calculations have been done for the three sources with 
different geometries. 

All the maps of IRAS 20I78-I-4046 consist of mainly 
one source, with a lobe seen in the ISOCAM maps. This 


source has been modelled with a single core and excellent 
fit has been obtained to the SED, radio flux and angular 
sizes. The dust density distribution which gives best fit to 
the data is of the form r“^ and the dust is dominated by 
silicates. 

While the IRAS HIRES and TIFR maps of IRAS 
20286-1-4105 consist of one source, ISOCAM maps show 
three sources at shorter wavelengths and two at longer 
wavelengths. This source has been modelled with two cores 
and excellent fit has been obtained to the SED along 
with consistency with radio observations. Dust in IRAS 
20286-1-4105 has been found to be coolest. The dust den¬ 
sity distribution which gives best fit to the data is of the 
form r“^ and the dust is dominated by silicates. 

The maps of IRAS 19181-1-1349 show multiple sources 
in all wavelengths. Hot spots have been found at positions 
away from the intensity peaks. Radiation transfer calcula¬ 
tions for this source have been done in cylindrical geome- 
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log [Wavelength (micron)] 


Fig. 14. (a) SED of IRAS 20178+4046. Symbols are - 
circles - flux densities from balloon-borne observations and 
IRAS HIRES, squares - from ISOCAM maps, triangles - 
from ground based observations of Faison et al. (1998) and 
dotted line - from IRAS LRS. The solid curve shows the 
SED corresponding to the best fit model, (b) Same as (a) 
but for IRAS 20286+4105. (c) Same as (a) but for IRAS 
19181+1349. The solid curve shows the best fit model for 
this source given in Karnik and Ghosh (1999). 

try with two cores. The dust distribution which gives best 
fit to the data is the one with uniform density i.e. r° and 
the dust is dominated by silicates. 

The differing density distributions tend to imply that 
IRAS 20178+4046, 20286+4105 and 19181+1349 are in 
progressive order of their evolutionary stage. 


Fluxes in four PAH bands have been obtained from 
the ISOCAM images. Ratios of the emission in different 
PAH bands have been obtained. Whereas two out of the 
three ratios have been found to be similar to the ratios 
obtained for other compact H ii regions the third ratio 
is quite different. Also, the band ratios for compact Hii 
regions have been found to be different from those for 
general Hii regions. Spatial distribution of PAH emission 
has been obtained and found to be diffuse and widespread. 
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Fig.11. The HIRES processed intensity maps of the region around IRAS 20286+4105 at (a) 12 /im, (b) 25 /rm, (c) 
60 and (d) 100 ^m. The contour levels are at 0.90, 0.70, 0.50, 0.30, 0.20, 0.10, 0.05 and 0.025 of the peak intensity. 
The peak intensities are 22.1, 175, 509 and 362 Jy arcmin“^ at 12, 25, 60 and 100 /rm respectively. 
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Fig. 13. Geometry of the models for radiation transfer calculations; (a) for IRAS 20178+4046, (b) for IRAS 
20286+4105. The dark circles in the centre represent embedded central stars and the grey areas the spherical shells 
of dust and gas. The white areas between the two contain only gas. For IRAS 20178+4046, the shell is heated from 
outside by the ISRF and from inside by the radiation from the embedded star. For IRAS 20286+4105, the radiation 
emerging out of the two inner gas-dust shells is combined and then used for heating the outer shell from inside. For 
more details refer to the text. 



